Tensorial representations of Reynolds-stress 
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The purpose of the present note is to contribute in clari- 
fying the relation between representation bases used in the 
closure for the redistribution (pressure-strain) tensor <|),;, 
and to construct representation bases whose elements have 
clear physical significance. The representation of different 
models in the same basis is essential for comparison pur- 
poses, and the definition of the basis by physically mean- 
ingfull tensors adds insight to our understanding of clo- 
sures. The rate-of-production tensor can be split into pro- 
duction by mean strain and production by mean rotation 
Pjj = P§.. +Pq u . The classic representation basis Q3[b,S,£2] 
of homogeneous turbulence [eg Ristorcelli J.R., Lumley J.L., 
Abid R.: J. Fluid Mech. 292 (1995) 111-152], constructed 
from the anisotropy b, the mean strain-rate S, and the mean 
rotation-rate D. tensors, is interpreted, in the present work, in 
terms of the relative contributions of the deviatoric tensors 



Pq ... Different alterna- 
5 (dev) 



tive equivalent representation bases, explicitly using Pi' 

'-> 

and Pq^. are discussed, and the projection rules between 
bases are caclulated, using a matrix-based systematic pro- 
cedure. An initial term-by-term a priori investigation of dif- 
ferent second-moment closures is undertaken. 



1 Introduction 

The pressure-strain (redistribution) correlation 



terms near solid boundaries |5|, obtained by the free-space 
Green's function solution of the incompressible flow Pois- 
son equation (linear in p') for the fluctuating pressure ||5). 
The most general approach to modelling in homogeneous 
turbulence is based on the order-4 tensors associated with 
the volume integrals (j)^' 21 ' and <|>y @, which require in 
particular, through obvious scaling arguments, that retained 

(r) 

models for (p.- should be linear in meanfiow-velocity gra- 

dients d Xj u~i EJ, and that <))•• should not depend directly on 
d^Ui 0. The redistribution tensor is symmetric by definition 
(JTj> and deviatoric, because of the incompressible fluctuating 

continuity equation = 2p'd X{ u'p = 0). Therefore second- 
moment closures (SMCs) represent redistribution as a linear 
combination of deviatoric symmetric tensors |4 6,[7). Let 



Sij :=| (d Xj Ui + d Xi uj) 
Clij :=j {d Xj u~j - d Xj Uj) i? d Xj u t - S t 
bij :=(u' t u' f ) u'lUj - |8y 
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Pij 

Djj := - pu\u'f) x M l ~ pu'ju' e d Xi U£ P s - P& (2g) 
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"plays a pivotal role in determining the structure of a wide 
variety of turbulent flows" flU. In (Q]), p is the pressure, 
Sij := ^(d Xj Uj + d Xj Uj) is the rate-of-strain tensor, are 
the velocity-components in the reference-frame with coor- 
dinates Xj, (•) denotes ensemble averaging, (•)' denotes tur- 
bulent fluctuations, f]).-;' 2 ^ is the rapid (meanflow/turbulence 



k (s;») 



is the 



interaction) volume-part of redistribution (4), (j)}. 
slow (turbulence/turbulence interaction) volume-part of re- 



distribution [4], and 



(„■) 



$ij + < l ) y stre me wall-echo 



viz the Reynolds-stress-anisotropy tensor bjj d2cb E), the 
production tensor Pij Oft appearing in Reynolds-stress trans- 
port [ID, which can be separated (Eft in production by mean 
strain-rate P$ (l2dT > and production by mean rotation-rate 



Pq, u (Eel . When homogeneous turbulence is submitted to ir- 
rotational strain Pq =0, while in case of solid-body rotation 



P§ ..= 0. The tensor Dy ( 2g i, which bears no particular name, 
appears in several early pressure-strain closures II2I8I9I . Naot 
et al. [8 1 attribute to Reynolds [2] its introduction in models 
for 0y. By ( Eg) , it follows that Dy can also be interpreted, 
indeed should be defined, as the difference of production by 
stain minus production by rotation. 



Almost all practical models^ for fl] El |3 gj 0121 E 
l9l [TTl[T2ll correspond to a linear combination of tensors in 
(fSjl and/or tensors constructed by products of tensors in (f2]) 
which are linear in the mean-velocity gradients. The par- 
ticular models range from simple phenomenological repre- 
sentations [ 1 , 118,2,9] using only tensors in (|2}, to more 
complex models making reference to irreducible representa- 
tion bases [4,6,7]. It is well established that the most general 
models can be represented in a basis of 8 linearly indepen- 
dent tensors [6J. However, there are several possible choices 
of the basis-elements, and the originally established forms of 
different models lfTTl2l[3Tl4Tl6ir7Tf8U9lfT11fT2l use different basis- 
elements. Although most of the models have been projected 
on a common basis [6|, the most complex one 1(131 1141 has 
invariably been expressed in a reducible form, including ten- 
sors which are linearly dependent [6] and can be projected 
on an 8-element basis lfT31l . To make detailed comparisons 
between different models, going beyond global evaluation of 
results against data for a given testcase, it is necessary to ex- 
press the models in a common basis. The purpose of the 
present work is to contribute in advancing towards the an- 
swer to the questions: 1) what is the common representation 
basis that should be retained, 2) what is the physical signifi- 
cance of the basis-elements, and 3) how different modelling 
choices associated with different routes followed in the con- 
struction of various models can be compared. 

In $2] we summarize results concerning the classical 8- 
element representation basis H, and debate on the argu- 
ments in favour of a polynomial representation basis vs a 
functional representation basis J4]|6]]. In <J3] we reinterpret the 
classical representation basis in terms of production by mean 
strain-rate P$ (ffdji and production by mean rotation-rate 
Pa-. d2el >. illustrating the physical significance to the 3 ba- 
sis elements, which represent all quasilinear models for rapid 
pressure-strain 1161 . In Sj4]we consider alternative bases built 
from the symmetric tensors [bij,Pij,Dij] or [bij,P§..,P[i..], 
which under the requirement of linearity in mean-velocity 
gradients are symmetric in [Pij,Dij] or [P^.j-Pq..], respec- 
tively, we show that the mean rate-of-strain tensor Sy can 
be explicitly projected on these bases, in which it does not 
appear explicitly, and discuss the advantages and drawbacks 
of such symmetric bases generated from tensors appearing 
in the transport equations for the Reynolds-stresses. In fJ5] 
we briefly note how projection of models on different bases 
can be made systematic using projection matrices and their 
inverses. Finally in SJ7J we illustrate, through a priori analy- 
sis of DNS data for fully developed plane channel flow, how 
term-by-term comparison of various models, expressed in a 
common basis, can be used to highlight different modelling 
strategies, indicating directions for future work (fjS]). 



'Notice that Launder 1 10 has suggested to further include the advection 
tensor Qj := pD(«J«^ in the representation, but this has not become standard 
practice. In the case of spatially evolving stationary quasi-homogeneous tur- 
bulence, Cij = piifd xi u'ju'j ^ 0, while in DNS studies of time-evolving spa- 
tially homogeneous turbulence, C/j = p3 f i(JHj ^ 0, implying that the sug- 
gestion of including C/j in the representation merits further study. 



2 Classical polynomial representation basis 

Recall IfTTl , that every order-2 tensor A £ E 3x3 
has 3 invariants, I A := trA, II A := 5(tr 2 A - trA 2 ) = 
\ ((A e( .) 2 -A em A hn ), and III A := detA = £^^2^3* = 
ZijkAuA ^Afs, coefficients of its characteristic polynomial 
K[jc] 9 p{x; A) := x 3 - l^x 2 + II A x - III A , whose roots 6 C 
(or G K if A is symmetric) are the eigenvalues of A, and 
which, by the Cayley-Hamilton theorem, is satisfied by A, ie 
A 3 = I A A 2 - II A A + III A I 3 El (4.24), p. 689]. Applying the 
Cayley-Hamilton theorem to (A + B) 3 - (A - B) 3 yields the 
Cayley-Hamilton theorem extension ABA = — A 2 B — B 2 A + 
IbA 2 + (I AB - I A I B )A - II A B + I A AB + I A BA + (I A 2 B - 
IaUb +IbII a )I 3 El (4.22), p. 688], and applying it to 
(A + B + C) 3 results in the extended Cayley-Hamilton the- 
orem ABC + CBA + BCA + ACB + CAB + BAC = (I B c - 
I B I C )A + (I CA - I C I A )B + (Iab - I A Ib)C + I A BC + I A CB + 
I B CA + I B AC + I C BA + I C AB + (I ABC + I CB a Wbc - 
IbIca+IcIab+IaIbIc)I3 El (4.21), p. 688]. The opera- 
tor Iq :— tr(-) denotes the trace (first invariant) of a tensor, 
and I3 E E 3x3 is the order-2 identity tensor in the Euclidean 
space E 3 . With respect to b d2cji, S (Eaji and d2bl . the 
Cayley-Hamilton theorem and its extensions [ 17] give ||4][6] 

b 3 = II b b + III b l3 (3a) 

bSb + b 2 S + Sb 2 = - II b S + I b §b + I b 2§ 1 3 (3b) 

bSb 2 + b 2 Sb=I b §b 2 + I b 2sb-III b S (3c) 

bfib= b 2 Q - Qb 2 - II b Q (3d) 

bfib 2 + b 2 fib= Illbfi (3e) 

The most general model for <|>y © is J3| a linear com- 
bination of b d2cj and b 2 + § Ufa 1 3 (I b 2 = — 2II b ) with coefi- 
cients which are functions of the invariants II b and III b , since 
by construction lb — (f2ct . The most general model for 
§ti^ O is S|6]|2) a linear combination of the symmetric 
deviatoric tensors constructed from products of btj (l2cl i. Sy 
(l2al ) and Ay ( |2bl , which form a representation basis II18II191 
of deviatoric-symmetric-tensor- valued isotropic functions of 
these 3 tensors, omitting, because of the linearity require- 
ment [3], elements nonlinear in the mean-velocity gradients 
(ie terms containing S"'J2" 2 with n\ +ti2 > 1). There are 
2 different approaches to constructing a basis: polynomial 
bases [18] and functional bases [ 19]. Polynomial bases IfTSl 
are formed by all products of integer powers of the generat- 
ing tensors which cannot be represented as a linear combi- 
nation of the basis-elements using CH-reduction, ie identities 
obtained from the aformentionned Cayley-Hamilton theorem 
and its extensions IfTTl , so that the representation coefficients 
are explicitly known polynomial (hence continuous) func- 
tions of the invariants. Functional bases [19| are potentially 
more compact, because they further reduce the elements of 
the corresponding polynomial basis by RE-reduction l20l pp. 
380-382], ie by solving appropriate linear systems [ 19 1. As a 
consequence, it can only be asserted lfT9l that the coefficients 
for representing a given product between integer powers of 



the tensors generating the basis as a combination of the basis- 
elements are functions (piecewise rational) of the invariants, 
not necessarily continuous 1191 . which are not always ex- 
plicitly known. For this reason |6| a polynomial representa- 
tion basis is preferable, and in the present case this is fT8ll6l 

«8[b,S,fi] := {f (B) , ne {!,■■■, 8}} with 



(2) 



(3) 



T 



(4) 



(5) 



(6) 



(?) 



(8) 



=b = fcyCj <8) ?j 

b 2 - = (btfbtj + |n b 8 )e/ 
=S := (ke _1 )%ei®e; 
zbS+Sb-Sl^b^-iPs-iS 

J bS 

=bfi - fib := (ke 1 ) (bfi - fib) ® ± P f 



=b 2 S + Sb 2 2 



I .* h 



3 V 



b 2 S 



=b 2 fi fib 2 
=b 2 fib bfib 2 



(4a) 
(4b) 
(4c) 
(4d) 

(4e) 
(4f) 

(4g) 
(4h) 



where the 8 tensors T/„) in (O were made nondimensional 
by scaling with appropriate powers of the turbulence kinetic 
energy k := ju' ( u' e and its dissipation-rate e := 2vd Xj u' j d Xj u' j . 

The above basis, <8[b,S,fi] := {t (n) , n £ {1, • • • ,8}} ©, is 
the classical [6| representation basis in homogeneous turbu- 
lence, where S and fi are usually fixed inputs of the prob- 
lem ED. 

Notice |6| that TVg) (l4hl can in principle be projected 
by RE-reduction to obtain an irreducible functional repre- 

sentation basis [0], Q3[b,S,fi] \ {T(g)}, but the coefficients 
for this projection are not known explicitly. Furthermore, 
even if the projection were to be sought, on a value-by-value 
basis |[20l pp. 380-382], the resulting representation co- 
efficients would not necessarily be continuous functions of 
the invariants |fT9l , and this makes functional representation 
bases awkward to use. The problem is even more acute in 
representation bases for wall turbulence [21|, where the ba- 
sic Gibson-Launder [22] rapid wall-echo model cannot be 

explicitly represented in the functional basis built by [b, S, fi] 
and the unit- vector in the normal-to-the-wall direction e„. 
Therefore, as noted in [6|, polynomial representation bases 
are the best choice, because the increased number of basis- 
elements is justified by the possibility of explicit continuous 
representation of models, using the Cayley-Hamilton theo- 
rem and its extensions ifTTl . 



3 Production by strain or rotation 

The tensor TV^ is related to strain-production d4dl i 

and the tensor T/^ to rotation-production d4eb . We 
can therefore construct an equivalent representation basis 
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Table 1. Matrix of coefficients Aht,,,,, for the representation H(„) = 
ELi fl HT„ m t (m) of the elements H (n) € Q3 (i ,)[b,S,P fi ;P s -] 
(|5]l as a linear combination of the elements T( m ) £ Q3[b,S,fi] 

\ 12 3 45678 
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Table 2. Matrix of coefficients a THnm for the representation T/ n \ = 
EJLi fl TH,„„H (m) © of the elements T (n) e <8[b,S,fi] © as a 
linear combination of the elements H( m ) £ Q3^/>)[b,S,p£j;P§] ©. 



|b,§,P fl ;Pj]:={H w ,n6{l,-,J 


I}} with 




H ( i).=b=. r (1) 




(5a) 


H (2) =b - jJ-b 2 ^ =• 1 (2) 




(5b) 


H( 3 ) -S =• T( 3 ) 




(5c) 


H (4) :=P s -:=(pe)- 1 (P % -|ft8 l7 ) 


e; ® ej 




|4c|4d| * 4 * 

- ^ 2 1 (4) - 3 1 (3) 




(5d) 


H (5) :=P fi :=(pe)- 1 P fl; / i ®^® 


2t (5) 


(5e) 



H 



(6) :=bP 5 + P 5 b-|l b * l 3 



(2l4fb}_4 1 -j- 



bS 

Q. + r fi 



S 

(i) +4II b T( 3 ) - |T( 4) +2T (6) 



H( 7 ) :=bPo + Pfib ^ 2T 



H 



(8) • = 



(7) 

2IIbt (5) 



2T 



(8) 



(5f) 

(5g) 
(5h) 



where the identities in (© are obtained by direct computation 
from definitions d2cH2ei >. using the following relations be- 
tween invariants, also obtained by direct computation from 
definitions $2^ 



L §1* St. B 

<4dj6bl i Py 

1 * 

bs 2 pe 

I8clgdl2i 15djgcl2> <4c|5dl2l 



(6a) 

(6b) 
(6c) 



l ^ ^ l ^ _ 4 j 4j ^ (6d) 

bD bP bP- b 2 S 3 bS 



bP 5 

(6Sl^_l T 

VS ~ 5 pe 4 bP s - 
EHU HTM) . 

b 2 P b 2 D b 2 P 5 , 

V9 3 ; pe 3 bP^ 



(6e) 



(6f) 



From the relations we can readily identify the matrix of 
coefficients a.. : = [ fl (HT)„,„] £ Kl 8x8 (Tab. Q]), and calculate 

its inverse a, . := a -1 (Tab. El, which relate the column- 

=(th) =ht 

* * * * 

vectors of basis-elements T := [T(ij, ■ ■ ■ ,T( 8 )] T and H := 
[H(i),---,H(8)] T © 



H = a T^T= fl -'H 

— =HT — — =HT — 



(7) 



a, . 

=(TH) 



The passage-matrix flj HT j (Tab. Q]), and its inverse a,, 
(Tab. EJ, can be used to systematically project models from 
one basis to another ($5). 



J,„:= P:= (pey^P^-jP^e^ej^ 

4T... IT... , OT._ ( 8c ) 



- f"(3) ~ 2 T(4) +2T( 5 ) 



j (4) :=D:= (pe)" 1 ^-^)?/®?/^ 5 
-fT(3)-2T(4)-2T( 5 ) 



(8d) 



J(5):=bP + Pb-|I rf l 3 

" 5^(1) +4n b f (3) - IT (4) +2T (6) + 2T (7) (8e) 
j (6) :^bD + Db-ll b6 . 3 ™ 

" f^d) +4n b f (3) - $T (4) +2T (6) -2T (7) (8f) 



j (7) :=b 2 P + Pb 2 -ll b2 .l 3 



21 _*T, 



VS J(1) + 5 I b | T (2)- 2III b T (3) 

- 2UbT(4) - 2IIbT( 5 ) - jT( 6 ) - 2T( 8 ) 



(8g) 



'(8) - : 



^D + Db 2 -|l b2 J 3 



mnu 



2II b T (4) +2II b T 



(5) " 3 ' (6) - 



2T 



(8) 



(8h) 



where the relations with the elements of 93[b,S,fl] := 

{T(„\, n 6 are obtained by direct compu- 

tation using the expressions © for the invariants, and the 
aformentionned Cayley-Hamilton theorem and its extensions 

ifTTll . Alternatively we can use P s .. I (p.. +Dy) and 

J2fl2 \ ^ 

Pq j(Pjj —Dfj) to obtain an equivalent representation 

basis generated by the 3 symmetric tensors b, P^ and P^, 
93[b,P s -,P fi ] := {F (n) , n € {!,••• ,8}} with 



4 Symmetric representation bases 

Many early models [2,8,9| used the rate-of -production 
tensor Pu (f2fb and the tensor D, -j which we identified in ( 2g 



as the difference between strain-production minus rotation- 
production. These 2 tensors are used by several researchers 
HUO, who sometimes [23 1 express models developed in 

the classical basis 93[b,S,fi] © in terms of P(j < f2fb and D, ; 
(2gi. If we consider the tensors P(j ( f2fb and D, ; (2gi repre- 
sentative of the influence of mean-velocity gradients d Xj Ui on 
the 0y ([TJl, since it is Pjj (Eft , and not Sy d2al > alone, which 
appears directly in the transport equations for the Reynolds- 
stresses 1H, we may construct a representation basis gener- 
ated by bij ( f2cb and these 2 tensors, 93 [b, P, D] := {-!(„), n G 
{!,-■•, 8}} with 



(2) 



IMS 



HE). 



(2) 



(8a) 
(8b) 



F(i) -=b =. T (1) 



(2b) . 



F (2) :=b 2 -Il b 2l 3 ^T 



(2) 



F (4) ; - p n 



ED 



H 



(4) 



H 



(5) 



(5):=bPj + Pjb-fl * l 3 ^H (6) 

s 

P3 



(6) ■= bP a + Pa b =• H (7) 
— k'b b k 2 

3 b 2 P f 



F (7) :=b 2 P s - + P s -b 2 -|l.^ I 3 



M 9T f ,2, f 

- 2III b T( 3 ) + 2II b T( 4 ) - |T( 6 ) 



"(8) 



b 2 P fi + P s b 



123 



H 



(8) 



(9a) 
(9b) 
(9c) 
(9d) 
(9e) 

(9f) 



(9g) 
(9h) 



where ( 9g i is again obtained by direct computation, using 
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Table 3. Matrix of coefficients a ITnm for the representation J(„) = 
Em=i fl JT„ m T (m ) ( fTOab of the elements J (n) € 03 [b, P, D] ® as a 

* * * — — 

linear combination of the elements T( m ) E 03[b,S,H] ©. 
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Table 4. Matrix of coefficients a FJmn for the representation F(„) = 
lLi fl FT„ m t (m) (DEB of the elements F (n) G 03[b,P s -,P A ] © 

* * * 

as a linear combination of the elements T( m \ 

Each of the bases 03[b,P,D] © or ©[b.Pj.Pjj] © is 
generated by 3 symmetric deviatoric tensors. Therefore, un- 
der the constraint of linearity in mean-velocity gradients ©, 

03[b,P,D] © is symmetric in P (I8cl2f1 > and 6 (18d|2g| >, and 

Q3[b,Ps,Pn] © is symmetric in P 5 (I5dl2dt and P fi (l5el2el 

On the contrary, the classical representation basis 03[b,S,fi] 
©, being built by the 2 symmetric tensors b ( f2cb and S ( f2ab . 
and the antisymmetric tensor fi (l2bl , contains the element 

T(8) := b 2 ftb — bfib 2 (PEEK but not the nondimensional de- 

viatoric projection of bSb +b Sb = I b §b + I b 2§b — IllbS, 
which is CH-reducible by ( f3cT >. As a consequence, the in- 
terpretation 05 (p) [b,S, P^; Pj] © of the classical representa- 
tion basis in terms of strain-production Pg, , d2dl i and rotation- 
production Pq ( f2eb is not symmetric in P^ (I5dl2db and P^ 
(I5il2el 

The basis-elements of03[b,P,D] © or 03 [b, Pj, Pq] © 
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Fig. 1. Locus, in the (Illb, —lib) -plane of the invariants of the 
anisotropy tensor 1 3 , of the condition — |llb = — Illb for which 

the representation Ti n \ = Lm=i fl TF„ m F( m ) dlObb of the mean 

strain-rate tensor S in 03[b,Pj,PQ] would be singular (the line 
— |llb = jj — Illb lies outside of Lumley's readability triangle 
I3ll24l . but includes the 1-C point). 

do not contain explicitly Sij d2ab . Hence, the practical use 

of the symmetric bases 03[b,P,D] © or 23[b,Pj,Pjj] ©, 
hinges upon their ability to represent models containing ex- 
plicitly Sij d2al i in their formulation, or, in terms of bases, 
on the possibility to project the classical representation ba- 

sis 03[b,S,fi] © onto these alternative bases. We may 
readily identify, from ( 14I5I6I8I ), the matrix of coefficients 
:= [fl( JT ) ] 6 R 8x8 (Tab. [3}, representing the column- 
vector J := [Jm,>> - ,J(8)] T ® in terms of the column- 
vector of the classical basis-elements T := \Jm\, ■ ■ ■ ,T(g)] T 
©, and from (14I5I6I9I ). the matrix of coefficients a. . := 

[a(ft) ] G ffi 8x8 (Tab. ©, representing the column-vector 
F := [F(i),-- - ,F( 8 )] T © in terms of the column-vector of 
the classical basis-elements T := \Jm\,--- ,T(g)] T ©, 

J = a T <^ t = a J (10a) 

— =JT — — =TJ— x ' 

F = a f T = a F (10b) 

— =FT — — =TF— x ' 

The inverse matrices a, , := a (Tab. H) and a, , := cT x 
=(tf) =ft v ■— ' =(tj) =jt 

(Tab. ©, which represent the column-vector of the elements 
of the classical basis in the 2 symmetric bases ( [Tol l, are read- 
ily obtained by straightforward inversion using symbolic cal- 
culation ||25l . The coefficients (Tabs. |5]6l l for the represen- 
tation of the elements of the classical basis 03[b,S,£i] © in 

03[b,P,D] © and 03[b,p5,P^] © are rational (hence con- 
tinuous except at poles) functions of the invariants. All of the 
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Table 5. Matrix of coefficients a T¥nm for the representation T(„) = Lm=i fl TF„,„F( m ) dlObb of the elements T(„) G Q3[b,S,fi] © as a 
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Table 6. Matrix of coefficients a Thlm for the representation T/ n \ = £ m=1 a T j,„„J(m) JlOab of the elements T(„) S Q3[b,S,fi] © as a 
linear combination of the elements J( m ) G *B[b,P,D] (|8). 



coefficients of (Tab. [5]l and (Tab. |6]l which could _|_ l ° k (fe-fe^ I -II b 8 



become singular have the same denominator ^ + |II b — Hit 



8 i 2tt TTT P £ 



1 



(Tabs. 06). It turns out that ^ + §II b - III b ^ 0, because + (ill b - §) — (Pj r - j^kSy) 

the line II b = 5 — |lHb lies outside of the realizability tri- ^ £ 
angle [3| of the invariants of b, II b and III b (Fig. Q3, with + 3-^(5),; + 2 F 0)U ( iia ) 



the unique exception of the 1-C point (II b = — j, III b = ^j). k _ p k 

Hence, the projection dTob is valid for any realizable |Lll (27 + 3^ b — g^v = ( — "TSpg — 4^ b p)^'J 

anisotropy tensor fey, except at the 1-C point. Using the rep- ^ 

resentation coefficients a( T F) 3 „, (Tab. [5]) and fl( T j) 3 „, (Tab. |6j, + 5 — (pnbij-\- ^IIb5j 7 ) 

the strain-rate tensor is explicitly represented in these bases ^ £ 

b y " " +(?nb-j)-(^-ifts y ) 

+(4-^-^(0,7-1^ 



^ + lH b - IIIb) ^Sy =( - ^ - \l *)fey i} 1* 

1 * 1 * 

4 J 0)ij + 4 J (&)ij 



1 1 

Ja)ii + zJ(S)ii (lit)) 



The representation ( II 1 al l of 5y in 58[b,Ps,P$}] © is more 

compact than the represenation d 1 lbb in Q3[b,P,D] ([8]), be- 

cause in the first case the strain-related terms (F«vy, F(5)ij, 

F(7)ij'< ©) are separated from the rotation-related terms 

(F(4)ij> F(6)ij< F(&)ij> ®)' wnereas m tne second case they 
are coupled because of the identities (I2f|2g|. 



5 Change-of-basis relations 

Since different model proposals in the literature use dif- 
ferent bases ll4ll6ll7ll8llMTinM3n4l it is necessary to estab- 
lish projection rules for the representation coefficients of var- 
ious models in different bases. Any expression (closure) for 
§ij can be represented equivalently in any of the bases, with 
linear relations between representation coefficients. Con- 
sider 2 bases with column-vectors of basis-elements A and 
B, related by 



— = =(ab)— 



B 



: a, , A 
=(ab)— 



(12a) 
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6 Isotropic limit 

Inspection of the matrices relating different bases 
(Tabs. [T3-|6j indicates that the model-coefficient d 1 3ab of 5;; 
changes depending on the basis used, and that in the sym- 
metric bases, <8[b,P,D] © and <8[b, P^, Pj}] ©, % can be 
represented as a linear combination of the elements of the 
basis (HB, so that the Rotta-CrowH constraint Il8ll9ll7l[l2l can 
be easily represented by a linear relation between the model 
coefficients c^ F ( 1 1 1 al 1 3bb or ( II lb!13bl i. The Rotta-Crow 
constraint 18 9.171 fl2ll requires that, at the limit of isotropic 
turbulence, the model for should recover the analytical 



solution lirrib-yo^/7 — ~m^ij — T?J v^y- F rom the represen- 
tation d 1 1 ab we readily have lim^o pk^, = —jP§... This 
suggests that an alternative expression of the Rotta-Crow 
constraint 0l9ll71[l2| is 



where the existence of the inverse matric a, ,:= a, 1 , £ 

=(ba) =(ab) 

R 8x8 is equivalent to the linear independence of the basis- 
elements, and hence to the fact that both A and B form rep- 
resentation bases. Denoting c^ A £ R 8 and c^ B £ R 8 the rep- 
resentation coefficients of <|> in each basis, we readily have 



(pe) > = £ ^ A „A (w) = cJ A A 



n=\ 



A g (AB) I=( g J UI) £^) T B (12b) 



lim 6; j 

b-+0 J 



fTT26l8l 



10 3 l J ~ 1A 



10 J s u 



(14) 



The practical conclusion from ( TPfl ) is that the Rotta-Crow 
constraint 18, 9,i7 HT2ll can be simply included in a closure 

expressed in *B[b,Pj,PQ] © as a constraint on the model- 
coefficient of P^ dev \ The implication of this analysis goes 



beyond the idea of implicit satisfaction 1 9 1 of the Rotta-Crow 
constraint, as in the isotropization-of-production model 0, 

* * 

suggesting that representation bases built using [P,D] or 
[Pg, Pq] need not contain explicitly Sij. 



proving that 



T 

: £(ab)-* 



(12c) 



matrices a, , £ 
=(ht) 

R 8x8 (Tab.|3), a 



FT) 

8x8 



ie that the representation coefficients in the 2 bases are re- 
lated by the transpose of the passage-matrix relating the 
basis-elements. We may therefore write, using the passage- 
~ 8x8 (Tab.[D, £ (TH) et 8x8 (Tab.©, g (jx) £ 
£ R 8x8 (Tab. |1, g (Tp) eI M (Tab.0, 
and §, . el 8xi (Tab. |5J, expressing the basis-elements 
( I4I5I8I9I ) of any of the bases as a linear combination of the 
basis-elements of another basis, 



(j),- 8 * 8 

-TZ = £ c ^J(n)ij = L C ^n H {n)ij 
P fc n=\ n=\ 



8 

I 

n=\ 



° 8 

C ¥n J (n)ij = L C ^n F (n)i. 
n=\ 



(13a) 



7 Discussion 

The basic closure for <|>y (slow return-to-isotropy, RI (TJ, 
and rapid isotropisation-of-production, IP [ 1 1 1) uses the ten- 
sors bjj (l2cb and Pjj < f2fb . which have a clear physical signif- 
icance. It was shown in the present work that the elements 

of the classic |6| representation basis Q3[b,S,Q] © can be 
represented in terms of production by strain P§, ( f2db or by 
rotation Pq.. ( f2eb . as can be interpreted models j8l|9]|7J us- 
ing their difference D[j = 1 



The widely known 

LRR [9 1 and SSG [4] models correspond to different weight- 
ings of the influence of P$,. d2dl i and P^ ( f2eb . the basic IP 
model [ 1 1 1 equally weights both, and the recent DY [ 12 1 pro- 
posal completely weights out Pq (l2el >. 

To highlight this remark, consider the a priori term-by- 
term analysis of different SMCs, using DNS data 11271 for fully 



2 Launder et al. | 9 | attribute to Crow [ 26 (3.6), p. 7] the behaviour of the 
rapid part of redistribution at the limit of isotropic turbulence, and this has 
been since repeated by several authors (71 1121 . Notice however Q5) that the 
same constraint had also been given by Rotta 1 1 p. 558]. 
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Fig. 2. A priori term-by-term analysis of SMCs (LRR [9], LRR-IP [9], SSG 0], DY [12]) for the components of the redistribution tensor 
eventually augmented by the anisotropy of dissipation (tyjj — 2pe/? £ , 7 := — p(e^ — |e)), with respect to DNS data [27] for the log 
and outer regions (y + > 30) of fully developed incompressible plane channel flow at friction-Reynolds-number Re Xw = 2003 [27]. 



developed incompressible plane channel flow, where the con- 
tribution of each basis-element to the closure of <|>y is consid- 
ered separately (Pig. [21). In this initial comparison, the mod- 
els were projected in the basis Q3(/>)[b,S,PQ;Ps] ([5]), which 
gives the most compact representation, because 3 of the 4 
models contain explicitly Sy I9,4,. 1T21 . Their comparison in 

the symmetric basis S^P^Pq] Q, will be the subject of 

future work. The contribution of bij =: ffrny to the models 

was separated to a slow part (pe ^Hrny; Fig. [2j and 

a rapid part (pe Cj 'Hn^; Fig. [2j> whose coefficient Cj 
is proportional to P^ ( |6bl >, ie dependent on mean-velocity gra- 
dients (the coefficient of this term is ^ in the SSG |4| and 
the DY [12 1 closures). For the weakly inhomogeneous (the 
strongly inhomogeneous near-wall, y + < 30, region was not 



plotted; Fig. |2|i pure shear flow studied, the four different 
SMCs (LRR (9), LRR-IP |9|, SSG f4), DY (ill) yield simi- 
lar global results, clearly indicating that they all include the 
anisotropy of dissipation (Fig. [2]). All of the closures use a 

(<b s;aJ ) * 

very similar coefficient for the slow bij term (pe cy Hmyj 
Fig.O, which is the basic RI model JTJ. Only the SSG |4| clo- 

sure contains a nonlinear slow term (pe cx; " (2)y)» which 
reduces redistribution by increasing (j) VI (in the present case 
(jixt < 0) and reducing § a (Fig.©. 

On the other hand, the different closures yield simi- 
lar global predictions for the rapid part of redistribution, 
but with different weights on each of the 4 basis-elements 
with rapid coefficients ^ (Fig. |2). The contributions of 

P^.. (pe cf V) H (5)ij ) and pf™ ] (pe cf^H^j) to the lrr 

and LRR-IP closures |9| are quite similar (Pig. O, the con- 



tribution of S v ^ (p£ 4* H(3)ij) to the LRR El c l° 



sure for (j)^ along with the contribution of H dev ^ yielding 
roughly the contribution of p| dev ' to the LRR-IP (9) closure, (ftji of (j); ,■ (21 



(to) 



of incompressible closures to compressible flows is obtained 
by replacing the mean strain-rate by its deviatoric projection 



Sjj — lSg£&ij everywhere, including in the definition 



for which ^ = 0. For the SSG ||4] closure the COntri- 
rt 1 " 93 ) * 1 

bution of Pq. (pe H^yj) is roughly ^ compared to 
the LRR and LRR-IP closures [9|, while the contribution of 
f| dev) (pe cf'^H^ij) is roughly the same (Fig. ©. The 
SSG |4| closure compensates the reduced contribution of Pq 

by the rapid contribution of bij (pe c, Hcuij), to yield a 
global prediction for the rapid part of redistribution similar 
to the LRR and LRR-IP closures ||9) (Fig. |2). Nonetheless 
this is a very different modelling choice, since the SSG JU 
coefficient of the rapid contribution of b,j is proportional to 



Pk 



ID 



■Ip_, and hence is related to a mechanism associated 



with strain-production P§.. and not rotation-production Pq.., 
implying that the SSG closure weights more the former than 
the latter, compared to the LRR and LRR-IP closures 0. 
Of course this is further amplified in the DY flT2l closure 
(Fig. |2]), which completely weights out the contribution of 

Pq (pe cij H( 5 jjj), further increasing the rapid contribu- 

tion of bij (pe c y { 'H {l)ij ). 

Establishing which is the best model is largely beyond 
the scope of the present note, requiring systematic compari- 
son for a wide variety of flows. The purpose of the present 
comparison is, on the contrary, to highlight the usefulness of 
element-by-element comparison of different models. 



8 Extensions 

To represent in strongly inhomogeneous cases, such 
as wall turbulence, let us assume, as do most models for near- 
wall turbulence 121], that we can identify (either geometri- 
cally or through gradients of local turbulence quantities) a 
unit vector pointing in the dominant direction of turbulence 
inhomogeneity, (|er|| = 1). By adding the deviatoric pro- 
jection of the tensor product (g) e^, viz f\ '■— ® eV, — II3, 
to the 3 basis-generators we obtain (under the condition of 
linearity in mean-velocity gradients and the fact that all pow- 
ers of T| can be represented as a linear combination of T| and 
I3, because by straightforward computation T| 2 = ^T] + |l3) 
14 more basis-elements [15|, which form, along with the 8 
elements of the quasi-homogeneous basis, a 22-element rep- 
resentation basis for <]);; in wall turbulence. As an initial ap- 
plication of these ideas, the 22-element basis obtained by 

combining T| with the 3 generators of *8(/>)[b,S,PQ;P<j] := 

|H(„), n £ {l,-- - ,8}} © has been used 1T31 to represent 
and compare all known families of single-point SMCs in a 
common basis. 

As usual 0[6) the extension of closures for to a 
reference frame rotating with angular velocity £2 RF is ob- 
tained by replacing everywhere £2y by the intrinsic (abso- 
lute) mean rotation-rate Wij :=Q i; — Ey-j^R^. The extension 



9 Conclusions 

By splitting the Reynolds-stress production-tensor Py 
( f2fb into strain-production P§. d2dl i and rotation-production 

" -|- Pq. ., the physical interpretation of the 



Pi j = 



tensor Dy ( 2g 1, present in the original expression of many 
closures both early [2,8, 9 1 and more recent Q[23), is ob- 
tained as the difference between the two mechanisms of pro- 

duction,Dy^Pj..-P Q ... 

* * 

The classical representation basis |6| <8[b,S,ft] © of 
closures for pressure-strain redistribution <|)y, can be replaced 
by equivalent representation bases, and model coefficients in 
different bases can be calculated in a systematic way, us- 
ing passage-matrices (fQ. Such representation bases, eg 

23[I>,P5,Pq] ©, are generated exclusively by tensors ap- 
pearing in the transport equations for the Reynolds-stresses, 
in which the rate-of-strain tensor Sy ( f2ab is weighted by 
the anisotropy tensor bij d2cb . In this basis, Sy can be ex- 
plicitly represented, using the Cayley-Hamilton theorem and 
its extensions flTl . as a linear combination of fey, P§.. and 
their powers/products which are linear in mean-velicity gra- 
dients (11 lab , and the Rotta-Crow constraint |T||26) appears 
as a limiting constraint on the model-coefficient of P$ . The 
practical implication of these results is that we can use the 

bases Q3[b, P^, P^] © or <B[b,P,D] © to compare differ- 
ent models, including those constructed in the classical ba- 
sis Q3[b,S,fi] ©, but also to construct general models for 
which do not contain explicitly the numerically stiff term 

Sij EH). 

An initial term-by-term decomposition on the basis- 
elements of different models (Fig.© illustrates how substan- 
tially different weighting of various basis-tensors can lead 
to very similar global results in weakly inhomogeneous pure 
shear flow. The detailed element-by-element analysis of dif- 
ferent models for various types of flows, using the bases 

«B[b,P s ,Pn] ® or Q5[b,P,D] ®, will be the subject of fu- 
ture work. 
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